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There is a well known mathematical similarity between two dimensional classical polarization 
optics and two level quantum systems, where the Poincare and Bloch spheres are identical math- 
ematical structures. This analogy implies classical degree of polarization and quantum purity are 
in fact the same quantity. We make extensive use of this analogy to analyze various measures of 
polarization for higher dimensions proposed in the literature, and in particular, the N = 3 case, 
illustrating interesting relationships that emerge as well the advantages of each measure. We also 
propose a possible new class of measures of entanglement based on purity of subsystems. 
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I. INTRODUCTION 

There exists a well known mathematical similarity 
between classical polarization optics in two dimensions 
and quantum two level systems. The Stokes vector and 
Poincare sphere on the one hand pQ are analogous to the 
Bloch vector and Bloch sphere on the other The 
Pancharatnam phase in classial optics systems [3] corre- 
sponds to the Berry phase in quantum systems 0]. 

Note that the point at the origin of the Poincare sphere 
represents a completely unpolarized beam of light, while 
any point on the surface of the sphere represents a com- 
pletely polarized beam. In the case of the Bloch sphere, 
the origin represents the maximally mixed state while 
any point on the surface represents a pure state. This 
suggests an additional analogy between the two systems, 
that classical polarization is analogous to quantum pu- 
rity, and measures of the two quantities should therefore 
be identical. This analogy between polarization and pu- 
rity has been discussed by some authors for a particular 
measure of polarization in higher dimension [5] . 

In this paper, we start by discussing in more detail 
the analogy between the classical and quantum cases, 
followed by a review of classical polarization in two di- 
mensions. We then move to higher dimensions, with 
particular attention to the three dimensional case. In 
doing so, we come across measures for quantum purity 
from quantum mechanics, namely the standard purity 
and von Neumann purity [5J. We also analyze measures 
of polarization in three dimensions due to Barakat [7], 
Friberg et. al. [8], and Wolf et al. [9]. We then proceed 
to compare these measures analytically, numerically and 
give them physical interpretations where possible. Our 
analysis adds to and clarifies much of the discussion on 
measures of higher dimensional polarization in the liter- 
ature [Mil]. 

Furthermore, we point out that the entanglement of 
a bipartite pure state can be thought of as the purity / 
polarization of a subsystem is traced out. This suggests 
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that using unconventional measures of purity, we can cre- 
ate new and interesting measures of entanglement. 



II. POLARIZATION OF BEAMS AND PURITY 
OF QUBITS 

A. Classical polarization states 

Consider a classical beam of light propagating in the z 
direction. Using the stochastic theory of electrodynam- 
ics, the magnitude of the electric field in the x and y 
direction are taken to be probabilistic ensembles given 
by Ei and E 2 respectively [HI [TTJ, [TH] . The polarization 
state of the beam of light at any point is given by the 
2x2 correlation matrix T^ 2 \ defined as follows, 

rg> = (E.E*), i = l,2 (1) 

where () denotes the ensemble average over different pos- 
sible realizations. We then define the polarization matrix 
$ as the normalized T, that is 

$ (2) =r (2) /Tr[r (2)^ (2) 

Therefore the polarization matrix <I>( 2 ) has unit trace by 
definition. In the following paragraphs, we suppress the 
dimensional superscript on the polarization matrix. Al- 
ternatively, the four normalized Stokes parameters, s^, 
can be used to represent the polarization state pQ. They 
are related to $ as follows, 

S , = Tr[$a»] = fyrf., (3) 

* = \s^, (4) 

where cr° is the identity matrix, and a , a 2 , and er 3 are 
the three Pauli matrices a z , a x , and a v respectively. Ein- 
stein summation notation has been used, i.e. repeated in- 
dices are summed over. Lowercase Latin letters run from 
1 to 2 (corresponding to the two Cartesian components 
of the field), while lowercase Greek letters run from to 
3. Since $ is of unit trace, sq is always unity. It can be 
shown that $ must be a positive matrix, which in turn 
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implies the following condition on the normalized Stokes 
parameters [20) . 



4 + s% + s\ < s 2 



1. 



(5) 



Immediately we see that eq.([9]) becomes the equation of 
a unit sphere, and taken as a position vector Sx,S2 and 
S3 describe a point on or within a unit sphere, called the 
Poincare sphere. The degree of polarization, P^ 2 ' of a 2- 
dimensional polarization matrix $ is derived by writing $ 
as the unique sum of two polarization matrices, one com- 
pletely unpolarized (i.e. a multiple of the identity), and 
one completely polarized [18] . The degree of polarization 
is then the ratio of "power" contained in the completely 
polarized matrix to the total power. It is given by, 



P {2) = a/1 -4det($). 



(6) 



Using eq.Q to write eq.Q in terms of the Stokes pa- 
rameters, we find that the degree of polarization is, 



P {2) = 



(7) 



That is, the degree of polarization is simply the distance 
from the origin of the Stokes vector within the normal- 
ized Poincare sphere. This makes intuitive sense and is a 
natural measurement, since the origin (where = 0) 
represents the completely unpolarized state, and the sur- 
face of the sphere (where P^ — 1) represents the set 
of completely polarized states, and the value of P for 
other states is "linear" in the distance metric within the 
Poincare sphere. 



B. Quantum two level system 

The most general quantum state is expressed in terms 
of the density matrix p, which contains all the statisti- 
cally observable information of the state. This matrix is 
positive, Hermitian, and of unit trace. In the case of a 
qubit, it is of dimension 2x2. We can write p as a linear 
combination of the identity and Pauli matrices as follows 



v 



r x Oi + r 2 er 2 + r 3 a 3 ), 



(8) 



where r = (ri,r2, f^) is the well-known Bloch vector [2J. 
Note that eq.pj is in fact identical to eq.Q. Moreover, 
since p is positive, we can also show that, 



rl < 1. 



(9) 



Therefore, the Bloch vector also lies within a unit sphere, 
known as the Bloch sphere. It is clear that in the 2- 
dimensional case, the quantum density matrix p is anal- 
ogous to the classical polarization matrix and that the 
Bloch sphere is analogous to the Poincare sphere, since 
they have an identical mathematical structure. We also 
see that the state at the origin of the Bloch sphere is 



the maximally mixed state whereas states on the sur- 
face of the sphere are pure states. Comparing this with 
the Poincare sphere, where the origin is a completely un- 
polarized state and the surface contains completely po- 
larized states, this suggests a direct analogue between 
quantum purity and and classical degree of polarization. 
Therefore, measures of the two quantities are for most 
purposes equivalent, and we keep this in mind in what 
follows. 



III. MEASURES OF PURITY AND 
POLARIZATION FOR N DIMENSIONS 

We proceed to introduce various measures of purity / 
polarization that have been suggested in the quantum 
mechanics and classical optics literature. Since purity is 
a property intrinsic to the density operator and invariant 
of the basis used, it should be invariant under unitary 
transformations. Therefore, one can always choose the 
basis where the density matrix is diagonal, and therefore, 
purity should be expressible as a function of the eigen- 
values of p alone, which we write as Ai > A 2 > ... > Xn, 
for an ./V dimensional system. 

If we write the purity as function of the eigenvalues, de- 
noted by II(Ai, Ajv), then we require that it be a real- 
valued function that is scaled such that it takes values be- 
tween and 1. It should take value 1 for a pure state and 
for the maximally mixed state; that is, 11(1, 0, 0) = 1 
and U(l/N, 1/N, 1/N) = 0, respectively. 



A. Standard purity 

In quantum information science, the common measure 
of purity for a quantum state p in an N dimensional sys- 
tem is given by Tr[p 2 ] [BJ. It takes a maximum value of 1 
for a pure state, and minimum value of for the maxi- 
mally mixed state. This purity is often scaled linearly so 
it varies between and 1, giving the following expression, 
which we call standard purity, 



n s (p)E 

In terms of eigenvalues, 



n s (/o) 



NTr[p 2 



1 



N - 1 



^Ef =1 A 2 -i 



N - 1 



(10) 



(11) 



B. Von Neumann purity 

Shannon entropy is used in classical systems to quan- 
tify uncertainty about a random variable. Von Neumann 
entropy generalizes this to quantum systems, and is given 
by, 

S(p) = -Tr[plog 2 (p)}. (12) 
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This measure quantifies the departure of a state from a 
pure state, i.e. its "mixedness" [6]. Note that the entropy 
of entanglement (a popular measure of entanglement for 
bipartite pure states) is defined to be the von Neumann 
entropy of one of the subsystems when the other subsys- 
tem is traced out. This implies that the von Neumann en- 
tropy is a good measure of mixedness. Therefore, one can 
define another measure of purity, Tl v (p) € [0,1], based on 
the von Neumann entropy, 



n»(p) = i 



Tr[plog 2 (p)] 
log 2 iV 



(13) 



Expressed as a function of eigenvalues, von Neumann 
purity is given by, 



n„(p) = i 



i 



N 



log 2 iV^ 



^Ailog 2 (Ai). 



(14) 



C. Polarization purity for N=2 

For the simple case of the two dimensional system, we 
have already seen that the classical degree of polarization 
is given by eq. ([6]) . Therefore, we have the two dimen- 
sional polarization purity given by p( 2 \ defined as, 



P^(p) = Vl-4det(p). 



(15) 



In terms of the two eigenvalues of p, this can be written 
as, 



P (2) (Ai,A 2 ) = y/i - 4AiA 2 
= Ai — A 2 , 



(16) 



where in the last equality, we used the fact that 1 = 
(Ai + A2) 2 . Moreover, as pointed out in the classical 
polarization case, P^ 2 ' is simply the length of Bloch vec- 
tor that represents p. That is, if we express the density 
matrix p as a Bloch vector as in eq.([8]), we find that po- 
larization purity simplifies to, 



P( 2 )(f) = ,/ r 2 +r 2 +r 2 = |^ 



(17) 



We wish to generalize polarization purity P^ to N > 2 
dimensions. Obviously, in principle there are an infinite 
number of ways to do this. However, only a handful of 
them have physical significance. In the following sub- 
sections, we discuss three possible generalizations that 
naturally follow from eqs. (15), ( fl6| , and (17) respec- 
tively. In particular, we pay particular attention to the 
N = 3 case, since it corresponds to classical polarization 
in three dimensions, a problem about which there is no 
consensus as of yet. 



D. Barakat heirarchy measures of purity 

In the case of an N x N density matrix p, Barakat 
has introduced a hierarchy of N — 1 purity measures [7] . 



These measures are defined by first writing out the char- 
acteristic polynomial equation of p as follows, 



det (p - 



AI) = \ N 
= 0. 



dA 



JV-l 



c 2 x N - 2 



+ (-i) N c n 

(18) 



The roots of this polynomial equation are clearly the 
eigenvalues of p (by definition). Each coefficient Ck is 
the sum of all possible unique products of k eigenvalues 
of p. That is, 



c k = y. nv 

l<U<...ifc<iV.7=l 



(19) 



For example, if iV = 3, then, 



C\ = Ai + A2 + A3 — 1, 
C'2 = AjA 2 + A1A3 + A2A3, 
C 3 = A1A2A3 = det (p). 



(20) 



In fact C\ — 1 and Cjy = dct(p) both hold for any N. 
Moreover, it can be shown that each Ck is expressible 
in terms of Tr(p m ), for m = 2,...fc, or alternatively in 
terms of the first k Casimir invariants of p (23] [24] . For 
example, for any N, we have, 

C 2 = [1 - Tr(p 2 )]/2, (21) 
C 3 = [l-3Tr(p 2 ) + 2Tr(p 3 )]/6. (22) 

Therefore, the Ck are invariant under change of coordi- 
nates. If p is a pure state (i.e. has rank 1), then all the 
Ck are zero, except for C\ which is always unity. If p is 
the maximally mixed state (all eigenvalues are j^), then 
Ck = ( 'fc ) j^ir where ( ^ ) is the binomial coefficient. 

With this in mind and noting that Ck coefficients 
themselves can be thought of as a measure of purity, 
Barakat then defines a hierarchy of measures of polar- 
ization given by B^ N \p) for k = 2, ...N. Requiring that 
P4 2) (p) collapse to P (2) (p) 

in eqJ15L one defines, 



B[ N \P)^ X 1- 



N k C k , k = 2,...,N, (23) 



The measure B^\p) takes the value zero for all k in 
the maximally mixed (i.e. the fully unpolarized) state, 
and takes the value 1 for all k when p is a pure (fully 
polarized) state. To get a feel for these measures, let us 
explore and simplify them using eq.(22) for some specific 
values of N and k. For N — 2, we have 



B ( 2 ] {P) = V / l-4det(p)=p( 2 )(p), 
as we required. For N = 3 we have, 



B?\p) 



3Tr(p 2 )-l]/2, 



Bf(p) = v / l-27det(p), 
= Vl - 27AiA 2 A 3 . 



(24) 



(25) 



(26) 
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For general N we have, 



F. SSKF purity 



B ( 2 N \P) 



NTi(p 2 ) - 1 



B 



N - 1 

N (p) = ^Jl-N"det(p), 

= y/l- N N X 1 \ 2 ...X N 



\W, (27) 



(28) 



Note the interesting relationship in eq.(|27[) where B, 
is simply the square root of the standard measure of pu- 
rity H s . However, B^ 1 ^ is unique, therefore we define 
Barakat's last measure of purity as lib, given by, 



U b (p) = B^\p) 



(29) 



We add lib to our collection of measures which will be 
compared to other measures later in this paper. 



E. KW purity 

Another measurement of 3-dimensional polarization 
proposed by Korotkova and Wolf relies on the total power 
in the completely polarized component [SI [TO]. That is, 
one splits the 3x3 polarization/density matrix into a 
unique positive linear combination of the identity ma- 
trix, a rank 2 matrix with degenerate eigenvalues, and 
a rank 1 matrix. The degree of polarization, Tlkw is de- 
fined to be the coefficient of the (fully polarized) rank 1 
matrix. We call this measure, the KW purity. It is given 
by, 



Hkw(p) = Ai - A 2 . 



(30) 



This is identical to eq.{ 16 1 in two dimensions. In fact, this 
particular measure has the same form for any N > 2, it 
is always the difference between the largest two eigen- 
values. The advantage of this measure is that it is very 
physically meaningful. It is the fraction of the power 
that can be rendered into a completely polarized field by 
passive linear elements. That is, if we use some (hypo- 
thetical) 3-dimensional polarizers, we should be able to 
extract the power of this fully polarized component. 

However, when one considers the rank 2 component of 
the polarization matrix (when broken up as we discussed 
two paragraphs above), one sees that this component is 
not fully polarized, but neither is it fully unpolarized. 
This suggests that it should have some intermediate po- 
larization of its own, and make a contribution to the level 
of polarization of the density matrix. Since ignores 
the rank 2 component completely, this suggests it is not 
well suited for as a measure of overall polarization, but 
is rather suited for measuring only the component of a 
field that is completely polarized. We clarify this in the 
next section with an illustrative example. 



There have been multiple suggested measures of three- 
dimensional polarization [STUOj. which rather than deal- 
ing with a beam propagating in one direction, deals with 
an arbitrary electric field distribution in 3 dimensions. 
The measure due to Setala, Shevchenko, Kaivola, and 
Friberg [5] starts by writing an equation similar to eq.p]), 
to decompose the 3x3 density matrix p. It is, 



1 



V3 



(31) 



where I is the identity, and Gi, i = 1 — 8 are a popular 
3 dimensional analogue of the Pauli matrices, known the 
Gell-Mann matrices [22], shown in the first Appendix. 
We have modified the coefficients from the Friberg paper 
to ease calculation. The eight coefficients r$ together 
form a (generalized) Bloch vector f. One can use eq.(31 ) 



together with the orthogonality and tracelessness of Gell- 
Mann matrices to show that, 



Tr[p 2 ] = 3 + 3|rf. 



(32) 



The density matrix condition Tr [p 2 ] < 1 implies that 

J2^=i r i — M 2 — 1) that is the Bloch vector r lies inside 
an 8-dimensional hypersphere of unit radius. The polar- 
ization is then defined in a manner analogous to eq.([7]) 
and eq.(17l as the length of the Bloch vector, i.e. the 



radial distance from the origin in this hypersphere, given 
by, 



H. 



*hf(p) = [ 



8 

»=i 



\f\. 



(33) 



We call this the SSKF purity or the radial purity, empha- 
sizing that it gives the length of a radial Bloch vector in 
a hypersphere. Equivalently, one can also invert eq. ( 32 ) 
to write the SSKF purity as a direct function of p, 



U sskf (p) = ^/[3Tr(p2)_i]/; 
or its eigenvalues Ai, A2, and A3, 



IWp) = \/~[3(A? + AI + A§)-1 



(34) 



(35) 



Note that the equation above is identical to eq. ( 25 ) , and 



therefore U ss kf{p) = B\ (p) for N = 3. For general N, 
we have, 



n sskf (p) = B { 2 N \ P ) = ^(pj : 

lNTr[p 2 } - 1 



N - 1 



(36) 



This is a very interesting property since each was osten- 
sibly derived in different manner. 
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The picture that has formed seems a natural gener- 
alization of the Poincare/Bloch sphere, since the centre 
of the 8-dimensional hypersphere still represents the to- 
tally unpolarized (or maximally mixed) state, and the 
states on the surface of the hypersphere are totally po- 
larized (pure). There is one essential difference however, 
that undermines the validity of this picture. In the 2- 
dimenstional case, all states in or on the Bloch sphere 
represent valid physical states and a positive density ma- 
trix. However, in dimensionality 3 or higher, it has been 
be shown that the physical constraint of positivity on the 
polarization matrix restricts the set of valid states to an 
irrelgular convex region that is a proper subset of the 
enclosing hypersphere. This physical region touches the 
surface of the enclosing hypersphere only in some places 
(where the fully polarized states lie) [231 HI]- That is, 
many states within the hypersphere and on its surface 
are unphysical since they would create polarization ma- 
trices that are not positive. 

The following diagram based on one by Kimura |23j 
shows all possible 2-dimensional cross sections of the 8- 
dimensional hypersphere, where the shaded regions rep- 
resent physical polarization/density matrices. Note that 
in fact most of the volume inside the hypersphere will 
be composed of disallowed unphysical states. A state 
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FIG. 1: Classes of cross section of the 8-dimensional space in 
which the generalized Bloch vectors live, based on diagram by 
Kimura [23]. In each diagram, the shaded region represents 
the allowable states, while the outer circle is the enclosing 
hypersphere. The pure states are where the shaded region 
touches the outer circle. Points A, B, C, D are specific states 
we examine. 



lying on the surface of the hypersphere is a necessary 
but insufficient condition for it to represent a fully po- 
larized physical state. States anywhere on the border 
of allowable (shaded) region must have at least one zero 
eigenvalue. A general state does not necessarily lie on a 
straight line between the maximally mixed state and a 
pure state. This must have been the case, since a pos- 
itive 3x3 matrix in general cannot just be written as 
a linear combination of the identity matrix (maximally 
mixed) and a rank 1 matrix (pure), there needs to be a 
rank 2 component as well. To illustrate these features, 
let us examine the states represented by points A, B, C 
and D in fig. ([I]). These four points are given by the 
following Bloch vectors: 

r A = (0,0,0,0,0,0,0,1), 
r B = (0,0,0,0,0,0,0,-1), 
r c = (0,0,0,0,0,0,0,1/2), 
f D = (0,0,^3/8,0,0,0,0,1/8). 



(37) 



Using eq.(31|, we can construct the corresponding den- 



sity matrices, and we find, 



PA 



PC 



2/3 

2/3 

-1/3 

1/2 

1/2 





PB 



PD = 









1/2 ' 
1/4 
1/4 



(38) 



We see that pa is not positive, and therefore unphysical, 
despite lying on the unit hypersphere. This illustrates 
the breakdown of the analogy with the two-dimensional 
Bloch sphere addressed above. The matrix ps however 
is positive therefore physical. Given that it is physical, 
we can see that it must be a pure state since it lies on the 
surface of hypersphere, and indeed it is. The matrix pc 
is also physical, and has a single zero eigenvalue, which is 
expected since it is at the boundary of allowable states. 
If it were to move slightly outside the boundary, the zero 
eigenvalue would become negative and therefore unphys- 
ical. The state given by pu is a typical state inside the 
allowable region. 

One may suggest that these properties are a result of 
an artificial asymmetry of the Gell-Mann matrices (in 
particular G3 and G$), and may be avoided if opt for 
a different basis set. However, this is not true, and the 
qualitative properties illustrated above are intrinsic to 
the N = 3 case, and still hold even if one exchanges the 
Gell-Mann matrices for a different basis set for SU(3) 
with the same basic properties of Hermiticity, traceless- 
nesss, and orthogonality. To see this, note that the sur- 
face of the hypersphere is 7-dimcnsional, and that pure 
states only form a 3-dimensional surface, we realize that, 
independent of the choice of the basis, the pure states 
form a very small part of the surface of the generalized 
Bloch hypersphere. Similar results hold for N > 3. 
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Yet despite the loss of the simple geometry of a filled 
hypersphere, the SSKF purity still, in some sense, quan- 
tifies the distance of the state from the maximally mixed 
state. Moreover, as we shall see later, it satisfies the 
intuitive depolarization criterion for this reason. 



IV. COMPARING PURITY MEASURES FOR 
THREE DIMENSIONS 

A. Graphical comparison 



Ai=0.2 




Thus far, we have discussed five contending measures 
of purity for N = 3 dimensions: the standard purity 
II S , the von Neumann purity Tl v , Barakat's last measure 
11;,, the KW purity Tlkw, and the SSKF purity measure 
Il ss kf- Since the standard purity II S is just the square of 
the SSKF purity H ss kf, we ignore the former and only 
include the latter in our comparison. 

To compare the four remaining measures, we set N = 
3, and recall that the purity will only be a function of the 
eigenvalues Ai, A2, and A3 = 1 — Ai — A2. Therefore, we 
plot the various measures of purity against A2 , for various 
set values of Ai. After examining the plots, we first note 
that rif, , H ss kf and n„ are always monotonically equiv- 
alent in the graphs. That is, they all increase together, 
decrease together, and have extrema at the same eigen- 
values. Put more specifically, derivatives of these three 
curves always have the same sign. In other words, for 
any two input states, these three measures always agree 
which state is more pure. 

If functions / and g are monotonically equivalent, we 
write / ~ g. Not that this is an equivalence relation. To 
rigorously prove that lib ~ H-sskf ~ n„ it suffices to show 
that each of them is monotonically equivalent to II S . It is 
clear that H ss k / = %/lTs ~ H s . In the appendix we prove 
that Ub ~ II S and n„ ~ II S . 

Returning to the graphs, we observe that the KW po- 
larization Tlkw often, but not always, gives the opposite 
results of all the other measures. This suggests that it is 
measuring something entirely different, and can be better 
understood through an example. 



B. Illustrative example 

Recall state C with spectrum {1/2, 1/2, 0} and state D 
with spectrum {1/2,1/4,1/4}. That is their respective 
density matrices are given by, 



Ai=0.33 



PC 



1/2 
1/2 




PD 



1/2 ' 
1/4 
1/4 



(39) 



Suppose we wish to use one of our measures of purity 
to find which density matrix pc or pu, is more pure. 
If we use the measure of purity II S or any of the other 
three in the same monotonic equivalence class, we find 
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FIG. 2: (Color online) Values of purity measures for vari- 
ous eigenvalues of a three dimensional density matrix. Each 
graph has fixed Ai, with A2 against the horizontal axis, and 
A3 = 1 — Ai — A2. In each graph, the upper solid curve in black 
is Barakat's last measure Hb, the upper dotted curve is the 
SSKF purity II ss fe/, the lower dotted curve is the von Neu- 
mann purity II„, and the thick solid curve is the KW purity 
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that state pc has higher purity than state D. If we use 
we find the opposite, state D is higher in purity. To 
find out which of these cases is qualitatively expected, 
we first note that in general, the more mixed a state is, 
the closer all the eigenvalues are to each other, with the 
extreme case being the maximally mixed state where all 
eigenvalues are identical, and the purer a state is the 
more a small number of eigenvalues should stand out. 

We can now reason that each of the density matri- 
ces pc and po have two identical eigenvalues (a "mixed" 
property), but in pc, the third distinct eigenvalue is fur- 
ther away from the identical two than in pp (that is, 
|0 - 1/2| > |l/2 - 1/4|), therefore p c should be more 
pure. Alternatively, we can reason that since both states 
have one eigenvalue of 1/2, then they are equal in this 
respect, and the other two eigenvalues should be the de- 
ciding factor in which state is more pure. The remaining 
eigenvalues for matrix pjj are 0.25, 0.25, these are iden- 
tical, and for matric pc are 0.5, 0, these are as different 
as possible. So we expect that pc must have higher over- 
all purity. Therefore any of the measures monotonically 
equivalent to II S , are well-suited for the general idea of 
overall purity. 

However, what if instead of overall purity, we are in- 
terested in the component of the density matrix that is 
completely polarized (i.e. the component that can pass 
through a three dimensional polarizer unchanged), then 
we see that, 



the following, 



PD 



1 





(40) 



That is, po has a nontrivial fully polarized component, 
the magnitude of which will be given by Hkw(pD) = 0.25. 
The matrix pc however cannot be decomposed in this 
way, has no fully polarized component that would pass 
through unchanged through a three dimensional polar- 
izer, and therefore Tlkw(pc) = 0. So we see that which 
measure of purity is more suitable depends on what we 
are interested in measuring. 

If we are interested in overall purity, we can further 
narrow down our choice of measure to the SSKF measure 
of purity Tl ss kf whenever our systems of interest involve 
depolarizing channels, a popular type of quantum noise 
channel. This is shown in the appendix. 



x — Rkw — Ai — A 2 , 
y=3(A 1 +A 2 -2/3) = l-3A 3 . 



(41) 



Physically, x is the KW purity, i.e. the fraction of the 
power that is in the fully polarized component, and y 
can be thought of as the fraction power that is not in 
the completely unpolarized component. In other words, 
x represents the power in the rank 1 component of the 
density matrix, while y represents power in the rank 1 
or rank 2 components, i.e. the power not in the rank 
3. Both x and y vary freely between and 1, with the 
condition that y > x (it follows from A 2 > A3). Then we 
we can express the eigenvalues in terms of x and y, 



Ai 



A, = 



A, = - - 



2 V 3 
2 V 3 

y 
3' 



■x), 
-x), 



(42) 
(43) 
(44) 



We can make use of these expressions to show that A 2 + 
A 2 + A 2 = (2 + 3a; 2 + y 2 )/6. Using this result along with 
cq.(35l, Tlsskf is expressible as, 



n 



sskf 



(45) 



This tells us that H ss kf includes the purity from x (i.e. 
Ilfeju) plus an additional component from y. Note that for 
a given x, the minimum value y can take is x, in which 
case Tlsskf = x = Tlk w . This is expected, since for y 
to equal x, this means there is no power in the rank 2 
component, and all the polarized power is in the rank 1 
component, so both measures agree. 



V. RELATION TO ENTAGLEMENT 
MEASURES 

The entanglement of a bipartite system is directly re- 
lated to the purity of a subsystem once the other sub- 
system has been traced out. For example, say we have 
a bipartite system of two qubits, A and B, given by the 
Bell state, 



C. Relationship between SSKF purity II ss fc/ and 
KW purity U kw 

We have already discussed the properties, strengths 



and weaknesses of II 



sskf 



and Tli 



It is of interest to 



find a simple relationship between them with aid of a 
pair of suitably defined variables. The following analysis 
is similar to results by Sheppard [14] . In the N = 3 
case there are only 2 degrees of freedom in setting the 
eigenvalues (since they must sum to unity). We define 



I*: 



|00) 



111) 



AB 



V2 



(46) 



This system is maximally entangled. If we define pa as 
the improper density matrix of the first qubit once the 
second one has been traced one, we get, 



p A = Ti B [|$)(*|] 

'1/2 
1/2 



(47) 



Note that pa is maximally mixed. If we had traced 
out system A and kept ps it would have been identi- 
cal. Moreover, if |4 r ) jlB where a separable state, then 
Pa would have been a pure rank 1 matrix. So, we see 
that maximal entanglement leads to maximal mixedness 
in the subsystem, and no entanglement leads to a pure 
subsystem state. This argument shows that the mixed- 
ness (one minus the purity) of a subsystem is a good 
measure of entanglement of the whole system. 

A common measure of entanglement for a bipartite 
system is entropy of entanglement E. It is defined as the 
von Neumann entropy S (i.e. a measure of mixedness) 
of a subsystem once the other has been partially traced 
out. That is, it can be written as, 



E(|*)) = i-n p& B (|*)<*|)]. 



(48) 



What if we replace IT in eq. ( 48 1 with another measure 
of purity, say Tl ss kf or lib? This would potentially give 
rise to another class of entanglement measures with dif- 
ferent properties, which could possibly be more relevant 
for some applications. 

For example we consider bipartite systems with N = 3, 
i.e. a system of two qutrits. Such a system has been 
studied by some authors, and even geometric descriptions 
developed to help visualize its entanglement |25j . Say we 
had the following two states, 



l*c) 



|00) 



111) 



V2 
= |00) 



1 ^ n >+H 



V2 

^|oo ) + |(M 



in) 



122) 



(49) 



(50) 



If we partially trace out one system in each of these two 
states, we are left with the subsystem matrices given by 
cq.(39) in the previous section. Then the we can apply 



the whole discussion regarding which measures of purity 
are more suitable, Tl s and its monotonic equivalents, or 
Hkw The question of which is more pure, pc or p^, 
can be asked differently, which is less entangled, \^c) ° r 

In this case, we can further appreciate why Hkw 
favoured po as more pure, because it favours \&c) as 
more entangled. Note that \^c) is a two dimensional 
Bell state in a three dimensional system, and one can 
think of it in a very specific sense as more entangled 
than \^d) since it is equal to a Bell state (even if it is of 
a lower dimension). One can also reason that \^c) should 
have higher entanglement since it has no separable com- 
ponent, whereas can be written as a combination 
of the separable state 1 00) , and the maximally entangled 
qutrit state. 

On the other hand, \*S?d) is more entangled than \^c) 
(and therefore pc purer than pd), since both have the 
same coefficient for the 1 00) state, yet \^d) has both 
a |11) and 1 22) component while \^>c) only has a |11) 



with no entanglement in the third level whatsoever. Once 
again, which state is more entangled depends on what 
exact properties we are looking for. 

All of this of course assumes the bipartite state |^)_4 B) 
is pure. In general, this state can itself be mixed and is 
written as a density matrix pab, which complicates the 
situation, and gives rise to the large and growing number 
of entanglement measures in the literature [26] . 



VI. CONCLUSION 



We showed that quantifying quantum purity for an N 
level system is equivalent mathematically to quantify- 
ing the degree of classical polarization in N dimensions. 
Then we described and analyzed multiple different mea- 
sures of purity, finding that which measure is more suit- 
able depends on what we are interested in measuring. 

In the more common case of measuring overall purity, 
we can use either the SSKF purity H ss kf, the von Neu- 
mann purity Tl v , Barakat's last measure of purity lib, or 
the standard purity IT S , all of which are monotonically 
equivalent. Of these, the measure most consistent with 
formalisms on quantum channels popular in the quantum 
computing literature is, H ss kf, given by eq.(33|. How- 



ever, all of these measures are monotonically equivalent, 
and any of them will suffice for simple comparison of 
purity between states. The standard purity Tl s in partic- 
ular is the simplest to use and the most common in the 
quantum mechanics literature. 

If instead we are interested only in the component that 
is fully polarized, then the KW purity Hkw is a more 
suitable measure, and will yield the strength of only the 
fully polarized part, discarding other components. It can 
be shown that for N = 3, H ss kf and are related in 
a simple manner once we add a variable to represent the 
second degree of freedom. 

Moreover, there is a direct relationship between the 
entanglement of a pure bipartite state and the purity of 
one of its subsystems once the other subsystem has been 
traced out. This can be used to give insight into measures 
of entanglement, and possibly create new entanglement 
measures based on various measures of purity. 



Appendix A: Gell-Mann matrices 

The Gell-Mann matrices (Gi,i = 1, 8) are a popular 
set of generators for the special unitary group SU (3) [52] . 
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They are given by, 
Gi 



1 

1 




G, = 



Gh = 



G 7 = 



1 
0-10 



-i 



1 

0" 

-i 

i 



G 2 = 



G 4 = 



G fi = 



G 8 = 



-i 

1 


1 



1 

o' 

1 

1 



1 

71 



1 
1 

0-2 



They are all Hermitian, traceless, and satisfy the or- 
thogonality relation Tr [GiGj] = 25ij, where Sij is the 
Kroncckcr delta. 



Appendix B: Monotonic equivalence of lib and II „ to 

i+ 

Assume we have an N dimensional state, with eigen- 
values Ai, Ajv, with the first N — 1 eigenvalues inde- 
pendent, and Xn satisfying, 



N-l 



X 



N 



(1) 



Then by eq.([TT]), eq.([T4]), and eq.(p8|, 

n^E^ + ^-E^) 2 ' ( 2 ) 

n _ 1 , EVogaA,- + (1-1^)1(^(1 
V ~ + lo g2 N 



U 2 b =l- A JV A 1 A 2 ...A JV _ 1 (1 - ]T X j)> 



(3) 
(4) 



where all sums over j in this section run from 1 to 
N — 1. We have squared lib since it does not affect 
monotic equivalence, and makes the calculation more 
tractable. The three expressions above monotonically 

equivalent if fj^S and have the same sign for 
each i — 1,2, N — 1. We have, 



SO, 

dX z 

dTL v 



an? 

dX t 



2A i -2(l-^A i ) 
2A, — 2Aat 

i _^[log 2 A i -log 2 (l-^A j 
1 

[log 2 Aj - log 2 Aj] 



(5) 



log 2 A 



(6) 



— — iV JV Ai...Aj_iAi + i...AAr_i(l — A? ~ 



A Ar Ai...A i _iA i+ i...Ajv-i(A i — Ajv). 



(7) 



We see that f^f 



and 



an 



-g^- both equal a positive number 
multipled by Ai — A at, and therefore will have the same 
sign. Also note that Ai — Ajv will always have the same 
sign as log 2 Ai — log 2 Xn, since Ai > A/v log 2 Ai > 

log 2 Ajv- Therefore the three measures ±I S , II^ and Ilf, 
are monotically equivalent, and will always give the same 
result if one uses them to find out which of two quantum 
states is more pure. 



Appendix C: Depolarizing channels as a criteria 



Consider a depolarizing channel, an important type of 
quantum noise [6|. It is a transformation which depo- 
larizes the input quantum state with probability p, (i.e. 
replaces it with I/A), and leaves it unchanged with prob- 
ability 1 — p. This channel is given by the following su- 
peroperator, 



S{ P )=p- + {l-p) P . 



(1) 



Applying the SSKF polarization measure of purity in 
eq.(36l to the output of this channel, we find, 



U ssk f(£(p)) = (1 -p)IL ssk f(p). 



(2) 



We observe that eq.Q has a very simple and intuitive 
form, showing that the purity simply scales down by a 
factor of (1 — p) after the state passes through the depo- 
larizing channel. This intuitive relationship does not hold 
for other measure of purity, even ones that are monoton- 
ically equivalent to Tl ss k / . This suggests lip is a measure 
with more physical meaning, and more relevant whenever 
depolarizing channels are in effect, as they often are. 
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